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Q ■ Abstract 

^ ■ Aspects of coherence and decoherence are analyzed within the optical Bloch equations. By 

■^ , rewriting the analytic solution in an alternate form, we are able to emphasize a number of unusual 

features: (a) despite the Markovian nature of the bath, coherence at long times can be retained; 

(b) the long-time asymptotic degree of coherence in the system is intertwined with the asymptotic 
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CN ■ difference in level populations; (c) the traditional population-relaxation and decoherence times, 

00 ■ 

^^ I Ti and T2, lose their meaning when the system is in the presence of an external field, and are 
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^^ ' replaced by more general overall timescales; (d) increasing the field strength, quantified by the Rabi 

frequency, O, increases the rate of decoherence rather than reducing it, as one might expect; and 



Q-i! (e) maximum asymptotic coherence is reached when the system parameters satisfy Q,'^ = l/(TiT2). 
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I. INTRODUCTION 

A quantum system whose dynamics is of interest is often part of, or coupled to, a second 
system whose dynamics is irrelevant. Examples include the internal quantum dynamics of 
a molecule in solution, a qubit imbedded in a solid, the translational motion of a single 
particle in a gas, the dynamics of one part of a molecule, etc. The overall dynamics of the 
total system is given by the Hamiltonian 

H = Hs + Hb + Vsb, (1) 

where Hs and Hb describe, respectively, the free evolution of the system of interest (that 
we henceforth refer to as the "system" ) and the extraneous degrees of freedom (termed the 
"bath"), and Vsb accounts for their interaction. Ideally, the system dynamics is obtained by 
solving the Schrodinger equation for the full Hamiltonian, H, and then averaging over the 
bath degrees of freedom. However, more often than not, this route is intractable. Hence, it 
is common to replace the full dynamics associated with Eq. (^ by an approximate master 
equation l|,y,|3| for the system density matrix, ps = Tr^ [p]- This type of equation provides 
a description of the two effects induced by the bath on the system: population changes and 
coherence loss. 

Recently, the degree of quantum coherence of a system has become increasingly impor- 
tant. For example, both the coherent control of molecular processes j^ |5| and quantum 
manipulations in quantum computing, quantum information and quantum cryptography 
lO; IZl IM Nl; rely upon the ability to keep the coherence of a system as well as to counter 
the decohering effects induced by the environment. One long-standing approach to reintro- 
duce coherence in a system that is interacting with a bath is to irradiate the system with a 
coherent electromagnetic field. 

In this paper we provide new insights into the coherence and decoherence in a 
paradigmatic two-level system interacting with a decohering environment and a resonant 
continuous-wave (CW) electromagnetic field. The model that we focus upon is the standard 
Bloch equation [Id] wherein the bath is Markovian, i.e., the coherence that is transferred 
to the bath is lost from the system forever. Using an alternate to the standard solution 
to this analytic problem [ll|, |l2|, we show that useful new insights emerge into the way 
in which the thermal bath and the external electromagnetic field interact to produce and 
sustain coherence in the system. Specifically, we emphasize that (a) despite the Markovian 



nature of the bath, coherence at long times can be retained; (b) the long-time asymptotic 
degree of coherence in the system is intertwined with the asymptotic difference in level pop- 
ulations; (c) the traditional population-relaxation and decoherence times, Ti and T2, lose 
their meaning when the system is in the presence of an external field, and are replaced by 
a more general overall timescale; (d) increasing the field strength, quantified by the Rabi 
frequency, Q, increases the rate of decoherence rather than reducing it, as one might ex- 
pect; and (e) maximum asymptotic coherence is reached when the system parameters satisfy 
n' = l/iT^T,). 

The organization of this work is as follows. A brief description of the two-level system 
is presented in Sec. IH] The asymptotic and time-dependent solutions to this model are 
discussed in Sec. 11111 Finally, the main conclusions derived from this work are summarized 
in Sec. ITVl 



II. THE TWO LEVEL SYSTEM 

The dynamical evolution of a two-level system influenced by a thermal bath can be 
modelled by means of the master equation 

dps i 



dt h 



H'{t),ps -nps, (2) 



where H'(t) = Hs + Hmt(t), and 7^ is a superoperator describing the evolution of the bath 
and its effects on the system (i.e., Hb + Hsb)- The evolution of the isolated (free) two-level 
system is determined by the Hamiltonian 



Hs= J]^.K)(« 



j=l,2 

and 



H,^, = -E{t) Y, dK)(j|, 



which accounts for the atom-field interaction within the dipole approximation, with E{t) 
£ cos{u!t + (f) and S being the strength of the electromagnetic field. 

The simplest model for dynamics of this type is the Markovian Bloch equation, 

dpij _ i 



dt h 



H'it),P -^Pr,r (3) 
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where 7^ is written in terms of the phenomenological relaxation times: Ti = Tj j and T2 = Tij 
{i 7^ j). In the absence of the electromagnetic field, Ti provides the timescale for changes 
in the system (eigenstate) populations, Pi^iit), with a rate given by Fi = l/Ti. Due to 
system-bath elastic collisions, there are random changes in the system phases that affect 
the off-diagonal terms Pij(t) {i 7^ j) and subsequently lead to system decoherence at a rate 
T2 = I/T2. 

In the most general approach the values of Ti and T2 are not constrained. For example, 



Skinner and coworkers have shown 
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IJ], using a non-Markovian model of a two-level 
system linearly and off-diagonally coupled to a harmonic quantum-mechanical bath, that 
T2 can actually be greater than 2Ti, with T2 = 2Ti in the weak coupling limit. In the case 
of the standard Bloch equation [lO|, however, one has 

2ri > T2 (4) 

in order to ensure that the reduced dynamics for the system always leads to completely 
positive maps of the density matrix [l5[, i.e., that Tr [p|] < 1 at any time. Thus, although 
the Bloch equation is mathematically well-defined for any values of Ti and T2, they lose 
their physical meaning when Eq. Q is not satisfied. Here, the condition ^ is retained 

throughout the study. ^ 

Equation can be solved within the rotating-wave approximation 16j by introducing 
the following change of variables Hf : 

i?i = 2Im[pi2] = -i{pi2 - P21), 

i?2 = 2Re[pi2] = P12 + P21, (5) 

-R3 = Pll — P22- 

Here, R3 is the difference in population between the two levels, and -Ri and R2 are the 
imaginary and real components of the off-diagonal density matrix elements. With these 



new variables, Eq. (0) can then be rewritten 18| in the standard form of the optical Bloch 
equations for a two-level system as 

dRi 



dt 
dR2 

~dr " 

dRs 
"df " 
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T2R1 + Ai?2 + ni?3, (6a) 

ARi - T2R2, (6b) 

fii?i-Fi(i?3-^3), (6c) 



where fl = 2di2S/h (with d^ = (l|d|2)) is the Rabi frequency with which the system 
oscillates between the two levels in the absence of a bath, and A = {E2 — Ei) / h — u is the 
detuning of the laser frequency 00 from the \Ei) —>■ \E2) transition. 

The quantity R3 in Eq. (jUcj) is the thermal equilibrium population difference to which R^ 
asymptotically relaxes in the absence of the external field, and is defined as 

^3 = i_^g-fe^2,i/fci3T' y'^) 

with U2A = {E2 — Ei)/h. R3 indicates the degree of mixedness of the reduced density matrix 
at temperature T and in the absence of the external field (i7 = 0). Two limits are therefore 
evident: (a) i?3 — > 1 if T — > 0, and (b) i?3 — > if T ^ cxd. 

III. RESULTS AND DISCUSSION 

A. Asymptotic Coherence 

Our interest focuses on the nature of the system coherence and its dependence on the 
system and field parameters. Two quantities serve as useful measures of coherence: the 
purity, X) and the interference contribution, (. The purity is given by 



X = Tr[p|] = pI-^ + pI^ + 2\pi2\ 
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= \ + l{Rl + Rl + Rl), (8) 

and depends on both level populations and interference. By contrast, the interference con- 
tribution, which we define as 

C = \Pu\' = \{Rl + Rl), (9) 

describes the coherence in the energy basis. Hence, though far from an ideal measure of 
decoherence, x is particularly useful due to its basis-independence. 

At long times, when the system reaches the equilibrium, all time derivatives in Eqs. (jHl) 
are zero, and the (asymptotic) value of the Ri becomes, for the on-resonance case (A = 0) 
emphasized below, 

Rf = 0, 

peq _ TiFa ~ 

-Kq — 77-:^ ;r^ -Ka- 



^■3 



T1T2 + n^ 



Consequently, inverting Eqs. (0) gives the reduced density matrix elements 



Q: 



eq _ 


/ eq^,* _ 1 


Taf 


^ 


12 
eq 


{P 
1 

2 
1 
2 


2i; 2 
1 1 


_l+TiT2fi2_ 
^3 1 


11 - 
eq 


1 ^^ 


5 


22 — 


/ l + TiT2n2_ 


) 



^i(2vP-7r/2)_^_^ 



(10a) 
(10b) 
(10c) 



where it is apparent that p^li ^ pS- Observe that 



^5^ = -^r2f]e^(2'^+-/2)[p^?-p^l] 



Pl2 



(11) 



giving a relationship between the coherence elements of the reduced density matrix and the 
population difference. 

Numerous features regarding the coherence at long times are evident from Eqs. (fTUIl . 
Obviously, the coherence is totally lost {( = 0) in the absence of the field (i7 = 0). Most sig- 
nificantly, as emphasized further below, the asymptotic extent of the coherence (as manifest 
in the value of P12) is directly proportional to -R3, the difference between the final populations 
of each state. Thus, for example, if the temperature is high, R^ = and C = regardless of 
the initial conditions. Similarly, at lower temperatures, and also regardless of initial condi- 
tions, the asymptotic value of ( can be nonzero, and the system can have therefore gained 
coherence due to the combined influence of the non-zero asymptotic population difference 
(related to the nature of the environment) and the field. 

Studies that focus solely on decoherence J8[ often neglect population-relaxation processes 
by setting Ti = 00. In this case [from Eq. (|10ap ]. all coherence is lost at equilibrium, with 
P12 — 0- Hence, long-time coherences that can exist in the case of finite Ti are missed. This 
reliance of long-time coherence emphasizes the significant interplay between the long-time 
population difference, i^g'^, the time that it takes to reach that limit (as manifest in Ti 
timescales), and the long-time coherence, p^2- 

Interestingly, coherence loss also occurs with large Q, as is evident in Eq. (jlOap . That is. 



Ip^I is an increasing function of ^2 until fi^ = l/(TiT2), at which point |p^2l = yT-ijTiRi/^- 
Increasing Vt beyond this value causes a decrease in the asymptotic coherence. Note that in 
the particular case where Ti is assumed very large, the asymptotic coherence is a decreasing 
function of Vt for almost all Vt values. 



B. Time Evolution 

The analytic solution to Eqs. (jH)) for the time evolution of the Ri is known [ll|. Here we 
rewrite it in a somewhat more enlightening form. 

Consider the case of on-resonance excitation (A = 0). Laplace transforming Eq. (jUc|) . 
and defining the transforms with over-bars, leads to the following linear system: 

{s + T2)Ri + nR3 = Rf\ (12a) 

(s + r2)i?2 = 4°\ (12b) 

(s + ri)i?3 - QR^ = rP + T1R3/S, (12c) 

where the R\ are the initial conditions for each of the variables. 
As seen from Eq. fll2b|l . the time-dependence of R2, 

R^it) = Rfe-^^\ (13) 

is readily obtained, since this variable is uncoupled from the other two. Therefore, for fixed 
T2 and Ti > T2, the dynamics is described by -Ri and -R3 beyond times on the order of T2. 
To obtain Ri and i?3, the value of -Ri resulting from Eq. ()12cp is substituted into ()12ap . 
giving 

[{s + ri)(s + T2) + n^] Rs = nR^ + {s + r2)i?f + ^"^^'^^' ^3, (i4) 

where the roots of the factor multiplying ^3 are s± = —a ± (3, with 



= £i±Ii v/(r2-ri)2-4^2 

a 2 ' ^ 2 

Depending on the magnitude of the discriminant of (3, three cases result. If |r2 — Fil 7^ 2^2, 
then 

i?3 = ^ + -^ + ^^, (15) 

S S — S_ S — S-i- 



with 



"» - iS* - rilTfP^- '"^' 



.2 d(0) 



slR^ + s„A + rir2i?3 



s_ s. 



slRf + s+A + rir2i?3 .... 

^' = .4^. - s.) ' ^^'^^ 



where A = nR[ + T2R3 + T1R3. The inverse Laplace transform of Eq. ()15p leads to 

R,{t) = D,+ {Bse-^' + Cse^') e""*. (17) 

Now, introducing Eq. (fTTjl into Eq. ()12a|l . we obtain for the imaginary part of the coherence 

Riit) = D,+ [B^e-f' + Cie"') e""*, (18) 

with 

TiVt ~ 

^' - -^r- ^'^ 

_ (£++£1) 

Therefore, for |r2 — Fil 7^ 21], both Ri and R2 decay exponentially. However, for |r2 — Fil < 
2fi, (3 is imaginary, and Ri and -R3 oscillate as they decay. Below, the time-dependent 
factors depending on (3 and that are responsible for the fine structure of the decays will be 
labelled Fi. 

Specific results for the case of "critical damping", |r2 — Fil = 2f2, can be obtained by 
considering the limit of the previous expressions when /5 — * (for which s± = sq = —a). 
In this case both i?i and R^ exponentially decay with rate a, although the F^ are linear 
functions of time. 

These results [Eqs. ()17|) and ()18p] allow direct consideration of the rates of falloff observed 
for the populations and coherences in the presence of a non-zero field. Specifically, if the 
field is strong {2VL > IF2 — Fi|) both the populations and the coherences decay with an overall 
rate a < F2 whenever Ti > T2, albeit with superposed oscillations. Thus, in principle, there 
is no distinction between the falloff rates of the diagonal and off-diagonal elements of the 
density matrix. This is also evident in %, where the long-time decay goes as 2a = Fi + F2, 
and it is not possible to clearly distinguish between two different timescales in its evolution. 
Nonetheless, since the Bi and C, depend on Q, the amplitude of the Fi can be greater than 
unity, and therefore may lead to slower x decay rates than 2a at early times despite being 
multiplied by the exponential. In the limit oi Q —>■ cxo, the Fi become highly oscillatory 
functions of time with decreasing amplitude about Fi = 1, and the decay rate exactly 
corresponds to 2a. This can be easily seen by substituting Eqs. (fT^. (fTTj) . and ffTH)) into 
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Eq. (jSJ. For example, for an initial state comprised of population in the state |1) (i.e., 
i?f) = Rf^ = 0, i?f = 1), one obtains 



_ 1 1 



2 



r2 — TiA _,__2 n^ , /r2 — Fi 



2(^^J sm^/?t+^^^)sin2/3t 



e-'"*, (19) 



which effectively approaches 



X = ^ + ^e-^-, (20) 



when i7 ^ oo (and therefore (3 -^ Q). 

If the field is weaker (2^2 < IF2 — Fi|) there are two system decay rates for Ri and R3, 
a±/9, with the smaller of the two dominating at longer time. In the parameter region where 
2Q < IF2 — Fi|, the decay rate increases with increasing ^2, evidently reaching its maximum, 
in this region, at critical damping. However, numerical evidence reported below shows that 
the decay rate continues its increase with increasing Q in the stronger field region, where 
2Q > IF2 — Fi|. Note that only in the limit of very small Q does the standard field-free 
interpretation {Ri and R2 decaying with rate F2, and -R3 with rate Fi) apply, since both B3 
and Ci approach zero in this limit. 

1. T = oo 

Consider first the case where T = 00. Here, the field is unable to beat out, at long 
time, the thermalizing effects of the bath, thus leading to zero asymptotic coherence. This 
happens even if the field is CW and is on for all time. The coherence dynamics on the way 
to the asymptotic result is of interest. As an indication of the non-intuitive nature of the 
associated relaxation timescales, consider the results for ( and x? shown in Figs, ^a) and 
(b). A grey scale is used to indicate the magnitude of ( and x as a function of t and fl 
for an initial state comprised of population in |1). A number of interesting features can be 
observed. First, as noted above, the decay rate increases with increasing Q, even beyond 
the case of 2Q = IF2 — Fi| ~ 1.33. Further, the onset of oscillatory falloff above this value 
is clearly visible. The (overall) falloff timescale, given by {2a)~^, is in this case ~ 0.375, 
making evident the inseparability of the meaning of Ti and T2 (i.e., one cannot observe two 
clearly separated decay timescales in x). 

The alternative perspective in Figs. Efa) and (b), where ( and x are shown as a function 
of log]^o^2 and time for fixed Ti = 2.5 and fi = 1 is also enlightening. In this case, the 



rate of decay of ( and x ^I'fi seen to fall off increasingly rapidly as T2 becomes larger, until 
|r2 — Til = 2i7, (here corresponding to logioT2 = —0.2). At this point, with increasing T2, 
the rate of falloff slows down and oscillates (these oscillations are better appreciated in the 
case of Q. The effect becomes increasingly exaggerated with larger Ti, as is evident from 
Fig. 121 that shows log ( for Ti = 10'*. Clearly, regarding T2 as the timescale for loss of 
coherence makes little sense in this context. 

Related insights, here into the reinterpretation of Ti, are provided by considering, for 
example, populations pn, as shown in Fig. HI Focusing on Figs. Ufa) and (b), we see that in 
the case where the field is weak (the solid curve for Q = 0.2) the rate of population relaxation 
does indeed decrease with increasing Ti. However, the value of Ti loses this qualitative 
meaning when the field strength is increased, as is evident by comparing the dashed curves 
and dotted curves, corresponding to stronger fields, in each of panels (a) and (b). That 
is, despite the fact that Ti differs in these two panels by a factor of 4000, the populations 
relax at qualitatively similar rates. Further, note that the increasing field strength tends 
to be associated with an increase in rate of the population relaxation. Similarly, Fig. EJ^c) 
shows that the rate of population relaxation depends on T2 as well as Ti. Notice that here, 
effectively, the fastest relaxation is achieved for r2 — 2Q, i.e., |r2 — Fil = 2Q [as one can also 
see in Fig. El^b)]. Clearly, the presence of the field imposes less meaning to the traditional 
values of Ti and T2. 

2. Finite T 

Additional significant effects arise when the temperature is not infinite and the asymptotic 
population difference, R^, is nonzero. Consider, for example, the extreme case were R^ = 1, 
shown in Figs. El^a) and (b), for ( and x^ respectively. Specifically, a comparison with 
Figs, m^a) and (b) shows large qualitative differences at all but the highest fl. In particular, 
two effects are evident. First, the nonzero asymptotic value of (, due to the nonzero R^, is 
clear. Second, significantly, both ( and x ^^re seen to show regions of Q where the function 
decays to the incoherent limit (zero in the case of (; one-half in the case of x)? but then 
reestablishes coherence to reach a long time value that contains coherence. This is clear, 
for example, in the case of ( for Q between approximately 2.1 and 5.4. Similar results are 
seen for x for a range of Q between 2.4 and 5.7. Hence, it would be misleading to suggest 
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that decay to the incoherent hmit iniphes no reestabhshment of coherence. For example, 
for Q = A, the system undergoes a coherence revival for a time interval At ^ 0.4, and for 
Q ~ 2.3 the revival takes place at t ~ 2.5 and the coherence remains permanently. Both of 
the above effects weaken with decreasing R^. 

Finally, Fig. IHl stresses the fact that the maximum amount of asymptotic coherence de- 
pends on the three parameters defining the system evolution. Thus, given Ti and T2, the 
value of the frequency for which one has maximum asymptotic coherence is Qr = 1/V^i^2 
[see Fig. Efb)], which leads to \p'^^\ = w p-^3/4. In particular, for the values used in Fig. 
Qr ^ 1.15. 

IV. SUMMARY 

Aspects of the ubiquitous Optical Bloch equations have been examined with a particular 
focus on the coherence of the system. The qualitative view that populations relax with rates 
of 1/Ti, and coherences relax with rates I/T2 is seen to be misleading when the system is 
irradiated by an external field with Rabi frequency Q. Similarly, increasing Q is shown to 
increase the rate of decoherence, contrary to simple intuition. Finally, the possibility that 
the coherence can recur to nonzero asymptotic values after having decayed to zero is noted. 
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FIG. 1: (a) ( and (b) x as a function of the Rabi frequency and time for Ti = 1.5, T2 = 0.5, 
and i?3 = 0. For all figures the darkest shading corresponds to 0.20 for ("-plots, and 20 increasing 
lighter shadings are used until ( = 0. For x^plots the darkest shading corresponds to 1 and 20 
increasing lighter shadings are used until x = 0.5. 
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time 




time 

FIG. 2: (a) Log^Q (" and (b) x as a function of log^g ^2 and time for Ti = 2.5, and R^ = 0. The 
upper limit of the ordinate scale corresponds to T2 = Ti/2, as in Eq. @. For all the logj^g C~plots 
the darkest shading corresponds to C, = 0.30, and 20 increasing lighter shadings are used until 
C = 0.0003. 
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FIG. 3: Logxo C as a function of log^g^a and time for Ti = 10^, and Rs = 0. The upper limit on 
the ordinate is as in Fig. |21 
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FIG. 4: Pii{t) for three different values of the Rabi frequency: fi = 0.2 (sohd hne), 0, = 1 
(dashed hne), and $7 = 5 (dotted hne). To compare, two different values of Ti are considered (with 
T2 = 0.5): (a) Ti = 2.5 and (b) Ti = 10^. In (c), pii{t) for three different values of T2 (with 
Ti = 10"^ and Q = 1): T2 = 0.05 (sohd line), T2 = 0.5 (dashed line), and T2 = 5 (dotted line). 
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FIG. 5: (a) C and (b) x as a function of the Rabi frequency and time for Ti = 1.5, T2 = 0.5, and 

i?3 = i. 
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FIG. 6: (a) %, (b) (, and (c) pn for three different values of the Rabi frequency: fi = 0.5 (solid 
line), r2 = r^r = 1.15 (dashed line), and 0, = 4 (dotted line). In all cases: R^ = 1, Ti = 1.5, 
Ta = 0.5. 
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